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Abstract 

The  influence  of  a  set  of  satellite  oscillators  on  the  response  behavior  of  a 
master  oscillator,  to  which  the  set  is  coupled,  is  of  fundamental  significance  to 
structural  acoustics  and  beyond.  The  focus  is  largely  on  the  induced  loss  factor 
that  the  satellite  oscillators  generate  in  the  impedance  of  the  master  oscillator. 
Much  of  the  research  work  performed  on  behalf  of  this  investigation  employed 
basically  sprung-masses  for  the  satellite  oscillators.  A  sprung-mass  is  a  primitive 
type  of  satellite  oscillator  and,  as  such,  limitations  are  imposed  on  the  range  of 
applicability  of  these  research  works.  In  this  paper  more  elaborate  satellite 
oscillators  are  introduced;  and,  especially,  a  wider  range  of  coupling  forms  and 
strengths  are  investigated.  A  number  of  new  insights  are,  thereby,  obtained.  In 
particular,  this  paper  is  to  facilitate  further  studies  of  the  relationships  among  the 
linear  impedance  analysis  (LIA),  the  energy  analysis  (EA)  and  the  statistical 
energy  analysis  (SEA).  These  studies  are  in  progress  and  are  to  be  reported 
subsequently. 
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I.  Introduction 

In  a  companion  paper,  designated  as  Reference  1,  the  complex  is  specified  in 
terms  of  the  impedance  Z®(<»)  of  the  isolated  master  oscillator,  the  impedance 
{imnij.) of  the  (r)th  isolated  satellite  “oscillator”  and  the  coupling  impedance 
{kr  I  ico)',  the  coupling  is  between  the  master  oscillator  and  the  satellite  oscillator. 
This  complex  is  sketched  in  Fig. la  [2-11].  In  this  complex,  the  satellite 
oscillators  remain  uncoupled  to  each  other.  The  purpose  in  the  present  paper  is 
to  introduce  a  significant  extension  in  scope.  In  this  extension  the  isolated 
master  oscillator  remains  the  same.  An  isolated  satellite  oscillator,  the  (r)  th,  is 
specified  by  the  impedance  Z^icai)  that  may  consist  of  both  a  mass  and  a 
stifBiess  control  term.  The  coupling  of  this  satellite  oscillator  to  the  master 
oscillator  is  specified  by  the  coupling  impedance  and  a  gyroscopic 

coupling  coefficient  (G^.)  [12,13].  The  coupling  impedance  Zcricoi)  consists  of  a 
mass  and  a  stiffiiess  control  term.  A  complex  of  this  type  is  sketched  in  Fig.  lb. 
The  mass  and  the  stiffiiess  control  terms  in  this  complex  are  related  in  the  forms 

K=K„(l+it1^)  ;  (K„IM)  =ml.  (la) 

kj.  —  (1  + / )  >  ^or  ^ —  (Oy  9  (1^) 

^cr  ^ocr  ^  ^cr  )  ?  ^^ocr  ^  ~  ^cr  ’ 


(Ic) 
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where  the  pairs  {M,K},  {mj.,kj.},  and  {m^  are,  respectively,  in  reference 

to  the  master  oscillator,  to  the  (r)th  satellite  oscillator  and  to  the  coupling 
between  them.  The  parameters  {tJo),  and  (7^^)  are  the  corresponding 
stiffiiess  control  loss  factors,  respectively.  To  complete  the  definition  of  the 
coupling,  the  mass  (Wc^)and  the  gyroscopic  coefficient  (G^)  need  to  be 
specified  [12,  13].  These  are  specified  through  the  coupling  parameters  which 
are  defined  in  the  normalized  forms 

^cr  =  Kr  /'"r)  5  gr  =  [G'r  /K  rn^)]  ,  (Id) 

respectively.  With  the  assistance  of  Fig.  lb,  the  linear  equations  of  motion  of  the 
master  oscillator  in  situ  and  of  a  typical  satellite  oscillator  in  situ  are  derived 
R 

^o(<^)  +  [^cr(p)  ~  (2a) 

1 

[Zr  (o))+Z^^(co)]V,(a))  +  [Z~  (6>)+G^  ]V^(a))  =  0  ,  (2b) 

respectively,  where  V^io))  and  V^(co)aie  the  responses  of  the  mass  (M)  of  the 
master  oscillator  and  the  mass  (/«^)of  the  (r)th  satellite  oscillator,  respectively, 
(R)  is  the  number  of  satellite  oscillators  that  are  coupled  to  the  master  oscillator, 
Ps((o)  is  the  drive  that  is  applied  externally  to  the  master  oscillator;  the  satellite 

oscillators  are  not  driven  externally,  Z^io))  is  the  impedance  of  the  master 


oscillator  in  isolation 
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Z°(fi))  =  (/®M)[l-W-^(l+i7o)];  y  =  (®/<»„)>  (3a) 

Z^(co)  is  the  impedance  of  the  (r)  th  satellite  oscillator  in  isolation 

Zj.{0))  =  {ia)mr)[l-iZrf  (1  +  iTJr  )]  ;  =  (®r  /®o)  5  >  y)’  W 

Zcr^co)  is  the  impedance  of  the  coupling  between  the  master  oscillator  and  the 
(r)  th  satellite  oscillator 

ZcrifO)  =  {iconic  ){m^r -i^crf  (l  +  ^7cr  )]  J 
f^cr  =  l^r)  \  ^cr  =  (<»cr  ^  ^cr  =  (^cr  /  j)  »  (3^) 

Z~j.{(o)  is  an  impedance  that  is  related  to  the  coupling  impedance 
namely 

^cr  (^)  =  ('  ^  )  [^cr  +  (^cr  f  ^  ^cr  )]  .  0^) 

and  (^crX  (^cr)’ (7cr)»  (G^)are  defined  in  Eqs.  (1)  and  (2)  [12].  The 
superscript  (o)  is  reserved  to  designate  quantities  that  pertain  either  to  the  master 
oscillator  in  isolation  or  to  satellite  oscillators  under  certain  and  definitive 
impositions.  Thus,  for  example,  from  Eqs.  (2a)  and  (3a),  one  may  state 

Z^,{co)V^{a))  =  P,{(D) ,  (4) 

where  V°{Q)y\s  the  response  of  the  master  oscillator  in  isolation.  By  a 
straightforward  algebraic  manipulation  of  Eq.  (2)  one  derives 


Z^(co)  V,(a?)  =  P,(a?)  ; 


V,(ci?)  =  B,(a>)V,(a?) , 


(5) 
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where 

R 

z„(0)  =  ZJ(®)  +  i;[{Z,(®)Z„(<!))}+(a,)2)  [Z/ffl)+  Z„(®)r' ,  (6a) 

1 

((d)  =  [Z~  ((d)  +  Gr  ][Z^  ((o)  +  ((d)T^  ,  (66) 

(.Qcrf  =  +  (G^)2  ,  (6c) 

and  the  quantities  Zg((D),  Zj.((o\  Z^^((o) and  Z~^((d) are  stated  explicitly  in  Eq. 
(3).  Indeed,  from  Eqs.  (1),  (3)  and  (6)  one  obtains 

Zo((d)  =  Zo(y)  =  0‘6>A/’)[l-(j)"^{[l-5(y)]+/[;7o+75(;')]}] ,  (7a) 

Br((D)  =  -[m^^  +  (z^^)^  (1  +i7Jcr)  -  i(gr  /j)]  • 

[(1  +  m^r)-(z,^)^  (l+ijj^^)Y^  ,  (7b) 

where 

(Xrr  f  (1  +  i  llrr  )  =  (^r  f  (1  +  i  Vcr  )  +  (^cr  f  (1  + « Vcr  )  >  (8a) 

R 

[sCz)  -  in.m  =  (y?  S  K  {[1  -  (zrf  (1  +  ‘nr)]  • 

1 

lm„  -(Zrr)^(l+incr)]-(ncr/y)^)- 

l(l  +  mrr)-(Zrrf(l+inrr)T'}  (8*) 
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and,  again 

Zrr  =  {Xrr  I  y)  \  Zr  =  I  y)  ’  ^cr  =  i^cr  f  y)  I 

i^cr  lyf  ^cr  (^cr  f  (1  +  «  7ct  )  +  iSr  >  yf  ’  ^cr  =  iOcr  )]•  (^c) 

One  notes  that  the  compound  coupling  parameter  (jq^)  is  a  ftmctional  of  the 
mass  and  gyroscopic  coupling  parameters  (mc;.)and  (f^),  respectively.  These 
coupling  parameters  are  defined  in  Eq.  (Id).  One  also  notes,  with  satisfaction, 
that  the  dependence  of  the  terms  in  the  sum  on  the  gyroscopic  coupling 
parameter  (g;.)  is  quadratic  so  that  the  sign  assigned  to  the  gyroscopic  coefficient 
{Gy)  plays  no  role  in  the  influence  of  the  individual  satellite  oscillators  on  the 
impedance  of  the  master  oscillator.  The  gyroscopic  coupling  is  in  quadrature  to 
both,  the  mass  and  the  stiffness  control  couplings. 

Examination  of  Eqs.  (7)  and  (8)  shows  that  the  normalized  impedance  that  the 
satellite  oscillators  collectively  induce  on  the  master  oscillator  may  be  cast  in 
term  of  the  two- vector  {5'(y),  (7?),  which  is  a  function  of  (R),  as 

indicated.  The  evaluation  of  this  two-vector,  however,  is  predicated  on  explicitly 
specifying  the  two-vector  {Xyy,  riyy}{R) ,  its  two  supplemental  components 
{Xy,riy){R)  and  {x^y,ri^y){R)znA  finally  assigning  the  compound  coupling 
ijqcr)-  The  two-vector  {x^y,  rjyy)  {R)is  designed,  for  the  sake  of  convenience,  to 
stay  fixed  with  respect  to  variations  in  the  index  (r) .  In  this  design,  the  springs, 
that  are  placed  on  either  side  of  the  mass  of  a  satellite  oscillator,  are  set  to  be 
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similar.  With  this  design,  Eq.  (8b)  is  evaluated.  In  this  paper  the  induced 
reactive  factor  5(:p)is  not  considered,  here  the  prime  interest  is  focused  on  the 
induced  loss  factor  %(y).  The  induced  loss  factor  rjsiy)  is  examined  for  a 
variety  of  coupling  forms  and  coupling  strengths  as  well  as  for  a  number  of 
values  of  the  modal  overlap  parameters  (Z>^)  and  associated  with  the  loss 
factors  77^  and  77^^,  respectively  [1].  By  and  large,  (Zv)is  set  equal  to  (^>c;.)and 
they  are  equated  to  (ft) ;  ft  =  ft^  =  .  The  simplifying  equalities  are  imposed 

without  considerable  loss  in  generality.  Using  Eq.  (8b),  the  exact  evaluations  of 
T],{y)  are  executed  for  three  values  of  (ft);  ft  =  (0.1),  (2.0)  and  (10) .  In  these 
evaluations  one  finds  that  77^(>')is  a  fimction  (b).  On  the  other  hand,  the  first 
order  approximation  of  rj^iy),  designated  by77/(>;),  that  is  derived  from  the 
replacement  of  the  summation  in  Eq.  (8b)  by  an  integration,  is  found  to  be 
independent  of  (b).  The  evaluations  are  graphically  displayed;  three  exact 
evaluations  for  the  three  values  of  (b)  and  the  corresponding  first  order 
approximation  of  ri^{y)dXQ  superimposed  in  each  of  the  displays.  The  com¬ 
parison  between  these  four  evaluations  can  thus  be  made  at  a  glance,  assisting 
greatly  in  the  interpretation  of  the  data  that  is  computed  and  displayed.  In  the 
frequency  range  of  concern,  when  (b)  is  small  compared  with  unity  the  levels  in 
the  exact  evaluations  of  77^  (y^)  undulate,  as  a  function  of  (y).  (The  frequency 
range  of  concern  spans  the  resonance  frequencies  of  the  satellite  oscillators.) 
The  excursions  in  the  undulations  increase  with  decrease  in  (b).  It  is  argued 
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that  the  first  order  approximation  constitutes  the  mean-value  averaging  of  file 
undulated  levels  of  rj^iy)  when  (b)  is  small  compared  with  unity  and,  again,  this 
mean-value  curve  is  independent  of  (b)  [1,  14].  As  (b)  approaches  and  increases 
above  unity  the  undulations  in  the  exact  levels  of  are  suppressed  and  the 

phenomenon  of  erosion  commences  and  increases.  In  the  fi'equency  range  of 
concern,  for  these  increasing  values  of  (b),  erosions  are  manifested  by 
progressive  decreases  of  the  levels,  in  the  exact  evaluations  of  rig{y),  fi'om  the 
levels  of  the  corresponding  first  order  approximation  [2,  3]. 
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11.  Resonance  Frequency  Distribution  of  and  Assignment  of 
Individual  Loss  Factors  to  the  Satellite  Oscillators 

It  has  been  argued  that  the  coupling  between  two  oscillators  may  cause  a  shift 
in  the  resonance  frequency  that  each  has  in  isolation.  Indeed,  it  was  suggested 
that  this  shift  may  be  used  to  determine  the  “coupling  strength”  [15,  16].  One 
may  then  question  as  to  what  exactly  are  these  shifts  and  are  they  significant. 
The  shifts,  as  such,  and  the  implications  that  they  may  harbor  are  not  addressed 
in  the  present  paper;  here  these  shifts  are  overridden  by  design.  The  design 
intends  to  derive  a  suitable  resonance  frequency  distribution  and  a  proper 
assignment  of  individual  loss  factors  to  the  satellite  oscillators.  The  design  is 
expressed,  then,  in  terms  of  the  two-vector  {x^^,  where  (x^^)is  the 

normalized  resonance  frequency  and  (7^;.)  is  the  loss  factor  associated  with  the 
(r)th  satellite  oscillator.  Examination  of  Eqs.  (7)  and  (8)  shows  that  the  nor¬ 
malized  resonance  frequencies  of  the  satellite  oscillators,  in  situ,  are  ascertained 
by  satisfying  the  equality 

(1  +  ~  (^/r)  >  0  ^c)  ~  (^rr)  ’  (^rr)^  ~  (^r)^ 

As  in  Reference  1,  here  too  (x^^)is  assigned  a  priori  with  equal  numbers  of 
resonance  frequencies  on  either  side  of  the  resonance  frequency  of  the 
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master  oscillator  and  the  distribution  is  aligned  in  ascending  order;  namely 

Xrr^Xqq',  q  =  {r  +  \)\  l<r<(i?-l).  {9b) 

A  scheme  that  simultaneously  satisfies  Eqs.  (9a)  and  (9b)  and  the  mid-point  just 
imposed,  demands  a  supplemental  condition  of  design  between  the  two  spring 
stififiiesses  that  support  the  mass  of  a  satellite  oscillator,  [cf.  Fig  lb.]  This 
condition  of  design  requires  that 

{Xj-rf  =  (1  +  fhcr)  for  7  =  (1/2) ,  (9c) 

where,  in  Eq.  (9c),  (F )  may  be  allowed  a  continuous  connotation  [1].  Consulting 
Reference  1  and  using  Eq.  (9a),  one  may  express  the  loss  factors  that  are 
associated  with  the  (r)  th  satellite  oscillator  in  the  forms 

rjr  =  [b^  I2{xrf]  dix^f  I  dr) ;  =  \b^j.  /2(x^^)^]  [5(x^^)2  Idr] ; 

Vrr  =  [(^rr)^r^  {(^r)^  Vr  +  i^crf  Vcr  (10) 

where  (/y)  and  (Z>^;.)are  the  modal  overlap  parameters  assigned  to  the  back 
spring  and  the  fore  spring.  These  two  springs  support  the  mass  (m^) .  [cf.  Fig. 
lb.]  For  the  sake  of  simplicity  and  computational  advantage,  the  stif&iess  on 
either  side  of  the  mass  of  a  satellite  oscillator  are  assumed  to  be  “similar”  in  the 
sense 

M  (^r  );  (^C/-)  =  (^r  )»  (^rr)  =  («r  +  (^r)>  (1 1) 

where  (a^)  and  (a^^)are  dubbed  the  spring  factors.  In  keeping  with  the 


definitions  of  (m^)  and  (g^),  as  the  mass  coupling  parameter  and  the  gyroscopic 
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coupling  parameter,  the  spring  factor  may  be  designated  the  stiffiiess 

coupling  parameter.  It  is  conducive  to  specify  (jc°)  in  a  form  that  is  compatible 
with  one  of  the  forms  introduced  in  Reference  1 ;  namely 

x“=[l+{l-2r}H*)^''^  (12) 

where 

r=/'(i?+i)~^  R  =  R{R+i)-^;  r(.R)  =  [rK2R)];  r<l,  (13) 
and  (r)  may  be  discrete;  1  <  r  <  ,  or  continuous  (s)<r<(R  +  e)  with  e<\. 
It  should  be  appreciated,  however,  that  although  in  this  paper  Eq.  (12)  is  cast  in 
stone,  other  forms  for  (x")may  be  readily  introduced  and  similarly  manipulated. 
Also  and  similarly,  as  introduced  in  Reference  1,  the  normalized  mass  (m^)of 

the  (r)th  satellite  oscillator  is  assumed  to  be  independent  of  (r)  and  to  be  of  the 
form 

R 

=  (m^  IM)  =  /M)(R)"* ;  M,  =  5](w^).  (14) 

1 

With  the  intended  exceptions  of  the  last  section  in  this  paper  it  is  convenient, 
without  a  great  loss  in  generality,  to  assume  that  the  spring  factors 
(«/•)  («cr)>  the  coupling  parameters  (f^)  and  (w^^)  and  the  modal  overlap 

parameters  (^y)  and  (^^^)  are  to  be  independent  of  (r) ;  namely 

Ctj.  =  a\  =  OCf.;  Sr  ~  ^cr  ~  »  W  ~  ^cr  ~  (15) 

and  it  is  observed,  as  already  intimated,  that  the  modal  overlap  parameters  (iy.) 
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and  are  set  equal  to  (b).  [cf.  Appendix  A,]  From  Eqs.  (9c),  (10),  (11)  and 


(12),  one  then  finds 

(«  +  a^)  =  (l  +  Wc),  (9d) 

Vr  =  Icr  =  Vrr=  T]{r)=(b  lit)  [K^)«)^  ]  ,  (16) 

where 

6=[(;rA)(5+l)-'].  (17) 

Under  this  imposition,  the  two-vector  {(x^^),  (rij^))  (i?)  assumes  the  simple  form 

(<7„)}  (R)  =  {(l+m,f^(x°),  !7(F)>,  (18a) 

and,  if  further,  the  mass  coupling  parameter  (fn^)  is  negligible;  (m^)  «  1 ,  then 

{(4)7  (7rr)}  (R)  =  {«)7  fj(r)h  im 


where  the  superscript  (o)  in  (x^^)  and  (77^^)  recognizes  that  (x®^)  and  (4)  arc 
restricted  to  specific  impositions,  [cf.  Eq.  (4).]  With  R=27  and  =  0,  Eq. 
(18)  is  evaluated  and  depicted  in  Fig.  2a,  as  a  function  of  (F) ;  Fig.  2a.  1  depicts 

(x^^)  and  Fig.  2a.2  depicts  In  Fig.  2a  the  modal  overlap  parameter  (b)  is 

increased  from  (0.1)  to  (2.0)  and  then  onto  (10).  To  these  changes  (7^;.)  increases 

by  a  factor  of  (20)  and  then  by  a  factor  of  (10  ),  whereas,  to  these  changes  in 
(Z>),  (x^^) remains  intact.  With  R  =  7andm^=0,  Eq.  (18)  is  evaluated  and 
depicted,  in  the  format  of  Fig.  2a,  in  Fig.  2b.  On  the  other  hand,  with  R  =  27, 
but  with  =0.75 ,  where,  for  example,  in  addition  a  =  1 .75,  =  0,  Eq. 
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(18)  is  evaluated  and  depicted,  in  the  format  of  Fig.  2a,  in  Fig.  2c.  It  is  noted  that 
for  a  continuous  r ,  except  for  obvious  end  conditions.  Figs,  2a.  1  and  2b.  1  are 
identical.  However,  Figs.  2a.2  and  2b.2  are  not  identical.  On  the  other  hand,  it  is 
noted  that  Figs.  2a.2  and  2c.2  are  identical,  however.  Figs.  2a.  1  and  2c.  1  are  not 
identical.  The  commonalties  and  disparities  among  Figs.  2a,  2b  and  2c,  are  as 
expected.  [It  is  also  noted,  in  passing,  that  and  (77^^) ,  stated  in  Eq.  (18)  are 
independent  of  the  gyroscopic  coupling  parameter  (g;.) .] 

It  remains  then  to  use  the  two-vector  stated  in  Eq.  (18)  to  evaluate  the  induced 
loss  factor  rfsiy).  Indeed,  using  Eq.  (18)  in  Eq.  (8b),  one  derives  the  more 
explicit  expression  for  in  the  form 


7isiy)  =  -{yflra 


R 


1 


-AmcUc  {z°f  {1+7 riir ) } -(g /  jp)^  }  • 


[il+m^)-{a+ac){z°f{\+iri{r)}}^Y  z°={x°ly),  (^d) 

where  (m^),  (g),  (y),  (a),  (a^),  (x^),  (m^)and  7(F)  are  stated  in  Eqs.  (Id),  (3a), 
(1 1),  (12),  (14)  and  (16),  respectively.  The  computations  of  ^^(y)  are  largely 
carried  out  assigning  the  standard  values 


(M,  I M)  =  10"^ ,  =  (0.1)  and  R  =  27  , 


(19) 
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where  (M^)  is  stated  in  Eq.  (14),  (b)  is  the  modal  overlap  parameter  and  (R)  is 
the  number  of  satellite  oscillators  in  the  set.  When  these  standard  assignments 
are  deviated  from,  specific  mentions  are  to  be  rendered,  notwithstanding  that,  at 
times,  the  employment  of  these  standard  values  may  be  reiterated. 
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ni.  Revisiting  the  Results  Presented  in  Reference  1 

It  may  be  useful,  at  this  stage,  to  reproduce  results  that  are  depicted  in 
Reference  1.  To  this  end,  the  following  impositions  are  rendered 

a=>0;  ac=l;  =  0;  gr=0.  (20) 

These  impositions  render  the  complex  commensurate  with  that  defined  in 
Reference  1  and  sketched  in  Fig.  la.  [cf  Fig.  lb.]  For  these  impositions,  Eqs. 
(16)  and  (18b)  and  Fig.  2a  are  validated.  One  finds  that  Figs.  2a.  1  and  2a.2  are 
akin  to  Figs.  3a  and  7a  of  Reference  1,  respectively.  In  addition  to  evaluating  the 
two-vector  specified  in  Eqs.  (16)  and  (18b),  the  corresponding  induced  loss 

factor  /;f(y)  is  evaluated  using  Eq.  (8d)  and  the  assignment  stated  in  Eqs.  (19) 

and  (20).  Again,  the  superscript  (o)  in  (7)  indicates  that  the  evaluation  is 

restricted  to  specific  impositions.  This  evaluation  of  7^  (;^)is  depicted  in  Fig.  3  a. 
The  modal  overlap  parameter  (b)  is  increased  fi-om  the  value  of  (0.1)  to  (2.0)  and 

then  onto  (10)  and  7/^(y)  is  evaluated  and  depicted  in  Figs.  3b  and  3c, 
respectively.  The  identity  of  Figs.  3a  and  3b  with  Figs.  5a  and  6a  of  Reference  1, 
respectively,  is  clear.  In  particular,  the  undulations  that  exist  in  Fig.  3a  and  the 
suppression  of  these  undulations  in  Fig.  3b  corresponds  to  a  phenomenon  that  is 
discussed  in  detail  in  Reference  1.  Figure  3c  does  not  have  a  counter  part  in 
Reference  1.  This  figure  is  included  in  order  to  bring  in  another  phenomenon; 
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the  phenomenon  of  erosion  that  may  beset  the  induced  loss  factor  (y)  [1-3].  It 

transpires  that  mean-value  averaging  of  levels  of  for  modal  overlap 

parameters  (b)  that  are  small  compared  with  unity,  coincide  and  are  thus 
independent  of  (b)  [1,  14].  This  coincidence  is  illustrated  in  Fig.  3d.  In  this 

figure  rjg  (y)  is  depicted  for  three  small  values  of  the  modal  overlap  parameter 
(6);  b  =  (0.01),  (0.1)  and  (0.3).  The  coincident  curve  representing  the  mean- 

values  of  T]g{y)  for  all  small  values  of  {b)\  b  « (1.0),  clearly  emerges  in  Fig. 
3d.  For  reasons  that  are  explained  subsequently,  the  curve  coincident  with  these 
mean-values  is  dubbed  (FOA).  An  erosion  is  here  a  phrase  to  describe  deviations 

of  the  levels  of  ri^{y)  from  these  mean-values  when  (b)  approaches  and 
increases  beyond  the  value  of  unity.  An  erosion  is  then  a  dependence  of  the 

levels  of  rj^  (y)  on  (b)  as  this  parameter  approaches  and  exceeds  unity.  To  bring 

into  focus  the  existence  and  nature  of  erosion  for  %  (y) ,  Figs.  3a,  3b  and  3c  are 
overlaid  in  Fig.  3e.  Also,  superimposed  on  Fig.  3e  is  the  curve  just  designated 
(FOA).  The  presence  of  erosion,  as  just  described,  is  thus  revealed  in  Fig.  3e.  It 
is  apparent  that  erosions  commence  at  the  edges  of  the  normalized  frequency 
range;  the  higher  the  value  of  the  modal  overlap  parameter  (b)  above  unity,  the 
more  the  inroad  from  the  edges  into  the  frequency  range.  (The  frequency  range 
of  reference  is  that  spanning  the  resonance  frequency  distribution  of  the  satellite 
oscillators.)  Figure  3e,  however,  shows  that  there  is  hardly  any  erosion  at  and  in 
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the  vicinity  of  y=l.  One  is  reminded  that  this  is  precisely  the  normalized 
frequency  region  where  (;;)  potentially  plays  the  more  significant  role  in 
controlling  the  response  behavior  of  the  master  oscillator.  To  cause  an  erosion  at 
and  in  the  vicinity  of  y-l,  the  normalized  overlap  parameter 

(Z));  b  =  [(;tZ>)(/?+1)“^  ],  needs  approach  and  exceeds  unity  [1-3].  This  extreme 
case  of  erosion  is  illustrated  in  Fig.  3f.  In  this  figure  the  only  change,  in 
parameters  that  specify  the  complex  assigned  to  Fig.  3e,  is  the  number  (R)  of 
satellite  oscillators.  The  number  (R)  is  (7)  instead  of  (27),  so  that  for  b  =  10, 

(b)  comfortably  exceeds  unity  in  Fig.  3f.  After  this  brief  digression,  it  is  time  to 
return  to  consideration  of  Reference  1  and  beyond. 

The  limited  scope  of  Reference  1  curtails  the  modeling  and  the  analysis  of  the 
complex  there  considered,  notwithstanding  that  complexes  employed  to  date  are 
largely  subjected  to  similar  limitations  [2-1 1].  A  poignant  question  arises:  what 
is  fresh  about  the  complex  defined  herein  as  compared  with  the  complex  defined 
in  Reference  1?  Whereas  in  Reference  1  a  satellite  oscillator  in  isolation  is 
characterized  by  a  mere  mass  control  term,  here  it  is  characterized  by  the 
oscillator  impedance  Z^(g>),  as  stated  in  Eq.  (3b).  This  oscillator,  in  addition  to 
the  mass  control  term,  may  also  possess  a  stiffiiess  control  term.  Moreover,  the 
coupling  impedance  as  stated  in  Eq.  (3c),  may,  in  addition  to  the 

stiffiiess  control  term,  possess  also  a  mass  control  term.  Finally,  the  coupling 
between  a  satellite  oscillator  and  a  master  oscillator  may  be  allowed  to  include  a 
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gyroscopic  control  term  [12,  13].  Obviously,  the  complex  sketched  in  Fig.  lb 
and  formulated  in  Eqs.  (5)  -  (8)  is  more  versatile  than  that  in  Reference  1,  as 
sketched  in  Fig.  la.  It  may  be  useful,  therefore,  to  investigate  a  few  of  the 
attributes  of  this  versatility  even  under  the  similarity  conditions  imposed  in  Eq. 
(1 1)  and  the  simplifying  assumptions  proposed  in  Eqs.  (14)  and  (15). 


Page  20 


IV.  Various  Coupling  Forms  and  Coupling  Strengths 

It  is  of  interest  to  evaluate  the  induced  loss  factor  77^(3;) ,  as  a  function  (y),  for 
a  variety  of  selected  coupling  forms  and  coupling  strengths.  The  coupling  forms 
are  defined  according  to  whether  the  coupling  is  dominated  by  either  a  stiffness 
coupling,  a  gyroscopic  coupling,  a  mass  coupling  or  combinations  of  these 
coupling  types.  The  coupling  strengths  of  these  forms  are  determined  by  the 
values  of  the  coupling  parameters;  the  stifBiess  coupling  parameter  (a^),  the 
gyroscopic  coupling  parameter  (g);  gr=g,  and  the  mass  coupling  parameter 
(m^) .  The  values  of  these  coupling  parameters  may  be  categorized  fi-om 
weaker-coupling,  to  moderate-coupling  onto  stronger-coupling  in  the  range  of 
values  of  0.03,  0.15  and  0.75,  respectively.  In  this  eategorization,  the  coupling  of 
the  satellite  oscillators  to  the  master  oscillator,  defined  in  Reference  1,  is  of 
stiffiiess  control  coupling  form;  i.e.,  g  =  =  0,  and  of  (very)  strong 

coupling  strength;  namely,  =  1.0  [a =0.0.].  Such  a  coupling  form  and  a 
coupling  strength  define  satellite  oscillators  commonly  designated  sprung-masses 
[1-11].  The  new  evaluations  in  this  paper  are  exhibited  in  Figs.  4-7.  Of 
significance  are  not  only  the  variations  in  the  coupling  forms  and  in  the  coupling 
strengths,  but,  also,  the  influence  that  changes  in  the  modal  overlap  parameters 
have  on  the  nature  of  the  induced  loss  factor  A  first  set  of  figures  is 

evaluated  with  6  =  (0.1),  a  second  with  b={2.0)  and  a  third  with  ^>=(10).  A 
major  feature,  common  to  all  evaluations,  is  that  the  undulations  in  the  first  set, 
for  which  Z>=(0.1),  is  suppressed  in  the  seeond  and  in  the  third,  for  which 
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Z>=(2,0)and  ^>=(10),  respectively  [1].  Figures  4-7  are  cast  in  the  format  of  Fig. 
3e;  corresponding  figures  in  the  three  sets  are  overlaid  so  that  the  undulations  in 
the  first  set  and  their  suppression  in  the  second  and  third  are  observed  at  a  glance. 
Another  major  feature,  common  to  all  evaluations,  is  that  the  erosions,  discussed 
briefly  in  Section  III  with  respect  to  1^s{y)  and  depicted  in  Figs.  3e  and  3f,  thread 
all  figures.  It  is  noted  that  when  b  =  [(;r6)(i?+I)"*]  approaches  and  exceeds 
unity,  erosions  occur  even  at  and  in  the  vicinity  of  These  extreme 

erosions,  however,  in  addition  to  the  dependence  on  (b),  seem  to  carry  a  slight 
dependence  on  the  coupling  forms  and  strengths.  These  dependencies,  which  are 
more  clearly  apparent  when  the  coupling  is  weak;  e.g.,  in  Figs.  6b  and  7d,  are,  at 
this  stage,  merely  noted.  Supplementally,  some  of  the  emerging  details  in  Figs. 
4-7  may  be  summarized  as  follows: 

1.  The  format  of  Fig.  3e  is  reproduced  in  Fig.  4a,  except  that  the  stiffhess 
coupling  parameter  is  reduced  firom  =1.0  to  =  0.75 ;  i.e.,  from  very  strong 
to  strong  stiffiiess  control  coupling  .  This  change  decreases  the  levels  in  Fig.  4a 
as  compared  with  those  in  Fig.  3e.  The  decrease  is,  however,  slight.  A  more 
drastic  decrease  in  levels  occurs  in  Fig.  4b  as  compared  with  Fig.  3e.  In  Fig.  4b 
the  stiffiiess  coupling  parameter  is  0.15,  whereas  in  Fig.  3e  ac=1.0.  This 
decrease  in  levels  is  largely  related  to  the  difference  in  the  coupling  strengths.  In 
Fig.  4b  the  coupling  strength  is  moderate. 
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2.  The  format  of  Figs.  4a  and  4b  is  reproduced  in  Figs.  5a  and  5b,  respectively, 
except  that  the  stifi&iess  coupling  form  is  changed  to  a  gyroscopic  coupling  form; 
namely,  =  ^^=0,  g  0.  In  Fig.  5a  the  coupling  is  strong;  g=0.75,  and  in 
Fig.  5b  the  coupling  is  moderate;  g=0.15 .  The  similarity  between  Figs.  4  and  5 
is  obvious.  Also  obvious  is  the  slope  in  the  curves  in  Fig.  5.  This  slope  is 
characteristic  of  the  gyroscopic  coupling.  The  gyroscopic  coupling  enters  in  the 
form  of  {gty)  and  not  merely  in  the  form  of  (g) .  [cf.  Eq.  (8).] 

3.  The  format  of  Figs.  4a  and  4b  is  reproduced  in  Figs.  6a  and  6b,  respectively, 

except  that  the  stiffiiess  coupling  form  is  changed  to  a  mass  coupling  form; 
namely,  a^=g=0,  .  In  Fig.  6a  the  coupling  is  moderate;  w^=0.15,  and 

in  Fig.  6b  the  coupling  is  weak;  fh^  =0.03 .  The  similarity  between  Figs.  4  and  6 
is  obvious  although  the  levels  in  the  former  are  higher  than  in  the  latter,  largely 
in  consequence  of  the  disparities  in  the  coupling  strengths. 

4.  Finally,  the  format  of  Fig.  4a  is  reproduced  in  Fig.  7,  except  that  the  stif&iess 

coupling  form  is  modified  to  accommodate,  in  addition,  another  form  of 
coupling.  In  Fig.  7a  the  additional  coupling  form  is  mass  control;  namely: 
ac  =  mc=0.75[a=0.1.],  g=0,  which  is  a  very  strong  coupling  strength.  In 
Figs.  7b,  7c  and  7d,  the  additional  coupling  is  gyroscopic  control;  namely,  in  Fig. 
7b:  =  0.53  [a:=0.43.],  g =0.54,  =  0,  which  is  a  strong  coupling  strength, 

in  Fig.  7c:  a^.  =  0.10  [a =0.9.],  g  =  0.1 1,  =  0,  which  is  a  moderate  coupling 

strength  and  finally,  in  Fig.  7d:  =  0.02[a=0.98.],  g  =  0.022,  ^^=0,  which 
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is  commensurate  with  a  weak  coupling  strength.  The  levels  in  these  figures  are 
set  largely  by  the  coupling  strengths;  the  higher  the  coupling  strength  the  higher 
the  levels. 

In  the  normalized  frequency  range  of  concern,  the  levels  in  Figs.  3-7  that 
pertain  to  modal  overlap  parameters  (b)  that  are  small  compared  with  unity;  e.g., 
b  =  (0. 1) ,  undulate.  On  the  other  hand,  again  in  the  normalized  frequency  range 
of  concern,  the  curves  in  Figs.  3-7  that  pertain  to  modal  overlap  parameters  that 
exceed  unity;  e.g.,  ^=(2.0)and  (10),  are  reasonably  smooth.  A  few  questions 
arise:  Can  these  features  be  estimated  by  replacing  the  summation  in  Eq.  (8d)  by 
an  integration  and  if  so,  can  this  integral  be  performed  with  ease?  Can  the  result 
of  this  performance  interpret  the  response  behavior  of  the  master  oscillator  in 
terms  of  the  various  parameters  that  define  the  complex?  And  last,  but  by  no 
means  least,  what  about  the  undulations,  do  they  feature  in  the  result  of  this 
integration? 
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V.  Replacing  a  Summation  by  an  Integration 

The  index  (r)  is  given  a  continuous  connotation  and  the  summation  in  Eq. 
(8b)  is  replaced  by  integration.  Under  a  condition  that  allows  this  replacement, 
Eq.  (8b)  assumes  the  form  [1] 

[i5(y)  -  irisiy)]  =  (y)^  ^ dz^(r)  U{r)n(j)]  • 

{[1  -  {z{r)f  (l+i?](r))]  [m,(r)-{z^iF)f  (1  +  -  [q,(F)/yf  }  . 

{[l+m^(F)]  -  [{z^iFyf  (1+/7^(F))]}"^  (21) 

where 

{qcir)lyf  =  Am,{F)  {z,(F)f  (l+i%(F))  +[{g{F)lyf- 
w(F)=//(F)(i?+l)-^  fiF)dz,{F)  =  dF',  {Zr{r)f={z(r)f+  {z^(F)}2; 

F  =  £-(/?+l)“^  ;  (^+F)  =  (j?+£')(7?+l)“^  (22) 

If  rjriF)  is  small  enough  so  that  the  vanishing  of  the  real  part  of  the  denominator 

in  the  integral  in  Eq.  (20)  predominates  the  values  of  the  integral,  the  integral 
yields 
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{S  (y)  -  i  T]{y)) ={7t  II)  /  [/(r^  )  nir^  )  ][1 + (''o  [{ ['"c  (''o  )  “  i^c  ('"o  )}^  ]  * 

{{z{ro)f  riifo)-  {Zc  {ro)f  tj^  )]+  4  (r  )  {z^  (7^  )}^  7c  (^  )} 

-i{fh,{r,)-{z,{r,)ff  +  [{z(r,)2,(^o)}"^(^o)^c(^c)]+[7c(^o)/3^f  }]  .  (23) 

where 

(^)<^  <(^+^); 

=  (1  +  ^cr)  =  {z(fo)f  +  Mro)f', 
7r„(^)=W^)}^7(^)  +  {Zc(ro)f  Vci^o) .  (24) 

Equation  (24)  defines  the  specific  value  of  (7^)  and  it  is  recognized  that  (7^)  is  a 
function  of  (y)  and  vice  versa.  Adopting  the  impositions  and  the  simplifications 
that  are  conducted  in  Eqs.  (16)  and  (18),  one  may  derive  from  Eqs.  (23)  and  (24) 
the  result 


7i,iy)=D[C  +  0{n(y)fV,  =  Z)C,  (25) 

with 

D  =  [^llr(RmM,  HMR)]  [(l+m<,)r' ,  (26a) 

C  =  mc+C!cf +(glyf'i^  (26b) 

<5  =  [(l  +  ^c-«c)«c]^ 


(26c) 
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and  where 

firo)={{a+acf'^^y{R){yfY^-,  lii{r^)={M,l MR),  (26d) 

vi.y)  =  Wo)  =  (b  l^)[Y{R){yf  ] ;  Wo)  =  7c(^)  =  %  (To),  (26e) 

(1+Wc)  =  (a+ttc) ;  [l+(r/2)r^^^^^  <  7  <  [l-(//2)r^^^2)  ^26f) 

The  underlying  condition  of  the  validity  on  Eq.  (25)  is  detailed  in  Reference  1. 
This  condition  holds  even  though  the  coupling  forms  are  elaborated  to  include 
not  only  a  weaker  stiffiiess  coupling,  but  also  mass  and  gyroscopic  coupling 
forms.  Moreover,  the  coupling  parameters  may  define  various  degrees  of 
coupling  strengths.  Strictly,  the  validity  of  Eq.  (25)  demands  the  equality  of 
^(j)^  the  primary  term  Tji(y)-  This  equality  designates  77/(;;)the  first  order 
approximation  to  ri,(y).  The  term  0{r]{y)f'^^  of  the  order  of  the  higher 
approximations  to  the  integral,  notwithstanding  that  situations  exist  in  which  (O) 
is  identically  zero;  e.g.,  when  the  satellite  oscillators  are  sprung-masses  for  which 
g  =  m^=0.  [cf.  Eq.  (20).]  In  these  situations  the  equality  of  the 
integral  evaluation  of  775(t^)to  rij{y)  need  not  be  specifically  invoked.  (More  on 
this  subject  when  in  a  subsequent  paper  higher  order  of  approximations,  than  the 
first,  are  to  be  evaluated.)  Clearly  and  significantly,  the  primary  term  rij{y)  and, 
therefore,  the  first  order  approximation  of  7j(y)is  independent  of  the  modal 
overlap  parameter  (b).  Without  much-a-do,  in  this  paper  the  first  order  approx- 
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imation  only  is  implemented  and  considered.  In  this  approximation  the  first 
order  approximation  to  rj^iy)  is  the  primary  term  7/(y) . 

Equation  (25)  reveals  the  parametric  composition  of  the  first  order 
approximation  to  the  induced  loss  factor  ri^{y)  =>  7/(y) .  Again,  one  is  remind¬ 
ed  tiiat  the  induced  loss  factor  (y)  describes  the  influence  of  the  coupled 
satellite  oscillators  on  the  loss  factor  in  the  impedance  of  the  master  oscillator. 
In  the  absence  of  couplings  this  loss  factor  is  {tJq)  ;  in  the  presence  of  couplings  it 
is  [^0+75(7)]-  Since  the  satellite  oscillators  add  merely  passive  elements  to  the 
complex,  rig{y)  is  invariably  positive.  Equation  (25)  confirms  this  statement  and 
exposes  the  proportionality  of  7/(y)to(D).  Therefore,  rjfiy)  is  directly  proport¬ 
ional  to  the  mass  ratio  (M^/Af);  it  is  estimated  that  for  a  reasonable  complex 
with  (Mj  /M)  equal  to  about  a  tenth,  (D)  is  of  the  order  of  unity,  [cf.  Eq.  (19).] 
The  quadratic  dependence  of  the  primary  term  rii{y)  of  rj^iy) ,  in  terms  of  (C) , 
which  entertain  the  term-components  (m^  +«c)  y)  5  is  of  significance. 

Again,  it  is  emphasized  that  the  primary  term  7/(y)is  the  true  first  order 
approximation  to  rj^^y) ;  in  this  context  the  term  DO{ri{y)}^  is  superfluous.  In 
any  case,  Eq.  (26)  indicates  that  this  term  in  ri^{y)  is  rarely  dominant  even  when 
the  loss  factor  ^(y)  of  a  typical  satellite  oscillator  exceeds  any  of  the  coupling 
parameters;  fhc ,  a^,  {g ! y) tjiy),  notwithstanding  that  in  the  absence  of  any 
couplings  rjgiy),  as  stated  in  Eq.  (25),  is  negligible  on  account  of  both,  C=0 
and  (9  =  0.  On  the  other  hand  when  couplings  do  exist,  in  both,  the  first  order 
approximation  and  the  exact  evaluations  of  %(y),  the  levels  are  certainly  not 
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negligible,  even  if  the  couplings  are  weak.  [Note  that  (7Jo)of  the  order  of 
(10“^)is  not  unreasonable  [12].] 

As  already  discussed  and  demonstrated  in  Figs.  3-7  the  curves  for  a  modal 
overlap  parameter  (b)  that  is  small  compared  with  unity;  b  «1,  possess  levels 
that  are  undulated.  The  excursions  in  the  undulations  are  the  more  pronounced 
the  smaller  is  the  value  of  (b)  [11].  [cf.  Fig.  3d.]  As  (^)  increases,  approaching 
and  exceeding  unity,  the  undulations,  again  true  to  form,  are  suppressed.  As 
already  intimated  in  Section  III,  increasing  (^>)  beyond  unity  brings  in  the 
phenomenon  of  erosion  which  worsens  as  the  modal  overlap  parameter  (b) 
reaches  higher  and  higher  above  unity.  A  remarkable  property  of  the  first  order 
approximation  (FOA)  of  [=7/(>^)]  emerges  when  this  quantity  is 

superimposed  on  the  respective  Figs.  4-7.  [cf  Figs.  3d,  3e  and  3f  ]  It  is  now 
observed,  in  these  figures,  that  the  mean-value  averaging  of  the  undulations  of 
the  exact  levels  of  7jg(y),  when  (b)  is  small  compared  with  unity;  Z>«1, 
converges  onto  the  first  order  approximation  (FOA)  of  Tj^(y)  [1,  14].  Since,  as 
already  observed,  the  first  order  approximation;  namely,  rjj(y),  is  independent 
of(Z>)  this  convergence  is  without  erosion  and,  of  course,  remains  so.  [cf.  Figs. 
3d  and  3e.]  The  independence  of  Tjj(y)  of  (b)  has  been  stretched  by  some  to 
conclude  that(^)  may  be  rendered,  a  priori,  equal  to  zero.  Neglecting  to  mention 
in  this  rendering  that  mean-value  averaged  levels  are  substituted  for  highly 
undulated  levels,  is  not  a  viable  scientific  procedure,  unless  ignorance  is  bliss 
[11].  On  the  other  hand,  when  (b)  approaches  and  exceeds  unity,  the  exact 
levels  become  free  of  undulations,  but  these  levels  erode  with  further  increases  of 
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(b).  The  erosion  in  the  exact  levels,  when  (b)  increases  beyond  unity, 
commences  and  progresses  from  the  levels  of  die  first  order  approximations.  To 
account  for  this  progressive  dependence  on  (b),  higher  order  approximations  are 
clearly  required.  What  is  doubtful  is  whether  higher  and  higher  order 
approximations  can  account  for  the  undulations  when  {b)  needs  to  remain  small 
compared  with  unity.  To  account  for  the  undulations  an  entirely  different 
approximation  procedure  is  thus  called  for.  In  this  paper  devising  such  an 
approximation  procedure  is  not  attempted. 
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VI .  A  Typical  Member  of  an  Ensemble  of  Complexes  Supporting 
Various  Parametric  Combinations 

In  the  preceding  evaluations  the  distribution  of  resonance  frequencies 
and  the  assigned  loss  factors  for  the  satellite  oscillators  are  sequen¬ 
tial  functions  of  the  normalized  index  (F) .  These  two  quantities,  exemplified  in 
Figs.  2a.  1,  2b.  1  and  2c.  1  and  in  Figs.  2a.2,  2b.2  and  2c.2,  respectively,  may  be 
smoothed  out  by  extrapolation  and  interpolation  into  monotonic  and  continuous 

functions  of  (F);  (F)=[r(i? +  1)~*  ].  This  kind  of  smoothness  is  rarely  found  in 
practice  and  a  question  arises  as  to  what  are  the  expected  consequences  of  more 
practical  assignments  for  these  parameters  and  others?  In  this  section  a  few 
layers  are  removed  in  the  quest  to  discover  the  phenomena  that  may  be 
encountered  in  the  induced  loss  factor  by  the  insertion  of  these  more 

realistic  parametric  values.  Since  the  assignment  for  the  parameters  that  define 
the  satellite  oscillators  and  their  couplings  to  the  master  oscillator  can  hardly  be 
drawn,  a  more  generalized  approach  is  undertaken  to  investigate  the  influence  of 
introducing  variations  in  these  parametric  values.  In  particular,  in  this  section 
two  parameters  are  selected  to  carry  these  variations;  either  individually  or  in 
unison.  In  the  first,  the  index  (r)  of  a  satellite  oscillator  is  assigned  a  pseudo- 
statistical  value.  [Pseudo-statistical  is  in  reference  to  a  sample  selected  out  of  an 
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ensemble  of  random  samples.]  The  index  (r)is  distributed  sequentially  and 
fractionally,  in  the  range  1  <  r  <  27.  A  pseudo-statistical  index  is  designated 
A(r),  where  A(r)  <  A(^);  ^  =  (r  +  1);  l<r  [cf.  Eq.  (9b).]  In  the 

second,  the  modal  overlap  parameter  (i’;.)is  assigned  a  pseudo-statistical  value 
that  is  distributed  in  the  ranges  (2)  >  (0.1)  and  (3.5)  ^b^>  (1), 

respectively.  The  distribution  of  A(r)  and  of  (b^),  with/?  =  27,  that  are 
employed  in  this  section  are  depicted  graphically  in  Figs.  8a,  8b,  and  8c.  The 
two-vector  {(x^^),  (7^^ )}(/?),  as  stated  in  Eqs.  (16)  and  (18),  is  typically 
depicted,  for  the  pseudo-statistical  values  shown  in  Figs.  8a,  8b,  and  8c,  in  Figs. 
9a,  9b,  and  9c,  respectively.  Figure  9a  depicts  (x^^)  as  a  function  of  A(r)  and 

Figs.  9b  and  9c  depict  (7;.^),  as  a  function  of  (F  ),  where  A(r)  =  [A(r)  (/?+!)“*] 

and  F=  [r(/?  +  l)“*].  [cf  Fig.  2a.]  It  is  observed,  in  Fig.  9a,  that  the  pseudo- 
statistical  variations  embody  the  phenomenon  of  mode  bunching  in  which 
variations  in  the  modal  density  of  the  satellite  oscillators  drastically  vary  as  a 
function  of  A(r)  [17].  On  the  other  hand,  as  Figs.  9b  and  9c  show,  the  loss 
factor  as  a  function  of  (r),  faithfully  follows  the  variations  assigned  to 

(bj.) .  In  Fig.  9b  some  of  die  values  of  (/y)  are  less  than  unity,  in  Fig.  9c  all  the 
values  of  (b^)  are  in  excess  of  unity. 

The  influence  of  the  variations,  described  in  Fig.  8,  on  the  induced  loss 
factor  ,  as  a  function  of  (y),  is  exemplified  in  Figs.  10  and  11.  Each  figure 
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represents  a  set  of  figures.  The  first  figure  in  each  set;  e.g..  Fig.  10a,  depicts  the 
base  situation  in  which  A(r)  =  r  and  =  1 .  The.  second  figure  in  each  set;  e.g., 

Fig.  10b,  depicts  the  situation  in  which  A(r)  is  as  shown  in  Fig.  8a  and  6^=1. 
The  third  figure  in  each  set;  e.g..  Fig.  10c,  depicts  the  situation  in  which 
A(r)  =  r  and  (^;.)is  as  shown  in  Fig.  8b.  The  fourth  figure  in  each  set;  e.g., 
Figure  lOd,  depicts  the  situation  in  which  A(7*)  =  r  and  {b^.)  is  as  shown  in  Fig. 
8c.  The  fifth  figure  in  each  set;  e.g.,  Fig.  lOe,  depicts  the  combined  situation  in 

which  A(r)  and  (b^)  are  as  shown  in  Figs.  8a  and  8b,  respectively.  Finally,  the 
sixth  figure  in  each  set;  e.g..  Fig.  lOf,  depicts  the  combined  situation  in  which 
A(r)  and  (bj.)  are  as  shown  in  Figs.  8a  and  8c,  respectively.  It  is  recognized 
then  that  each  set  of  figures  presents  a  complete  evolution  in  the  process  of 
applying  the  pseudo-statistical  variations  depicted  in  Fig.  8  to  the  two  parameters 

A(r)  and  (Zy).  Also,  each  set  selects  a  specific  coupling  form  and  a  specific 
coupling  strength.  Thus,  Fig.  10  pertains  to  a  strong  stiffiiess  coupling: 
a^=l.O  [a =0.0.],  -g  =  0,  and  Fig.  11  pertains  to  a  mix  of  stiffiiess  and  of 
gyroscopic  coupling  of  moderate  strength:  =  0.1[a  =  O.9.],  g  =  0.1 1,  =  0. 
[cf  Figs.  3e  and  7c.] 

The  first  figure  of  each  set;  namely,  Figs.  10a  and  11a  exhibit  undulations 
in  the  levels  of  the  induced  loss  factor  ri^{y),  as  a  function  of  (y).  However, 
these  undulations  are  small  and  they  are  completely  suppressed  as  soon 
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as  (6^)  approaches  the  value  of  (2).  [cf.  Fig.  3b.]  The  variations  depicted  in  Fig. 
8  are  clearly  discernible  in  all  the  subsequent  figures  in  the  series  entitled  Figs. 
10  and  11.  True  to  form,  there  is  but  a  tinge  of  edge  erosion  in  Figs.  10a  and 
11a,  yet  in  the  likes  of  Fig.  3b  stronger  sign  of  erosion  has  already  reared  its 
head.  From  Figs  10a  and  11a  to  Fig.  3b,  (b)  is  changed  from  unity  to  merely 
two.  To  confirm  this  statement  and  to  provide  for  convenient  and  interpretable 
data  from  which  to  judge  the  more  erratic  data  that  incorporates  the  pseudo- 
statistical  variations,  the  first  order  approximation  of  ,  given  in  Eq.  (25),  is 
superimposed  on  Figs.  10a  and  11a  and  on  all  other  figures  in  the  series  entitled 
Figs.  10  and  11. 

Again,  the  pseudo-statistical  variations  are  defined  by  two  competing  and 
nearly  independent  factors,  i.e.,  by  A{r)  and  ^(a')[=(^v)].  In  Figs.  9a,  10b  and 
11b;  it  is  observed  that  at  a  mode  bunching  (a  rich  modal  density)  region  the 
influence  of  the  satellite  oscillators  is  more  pronounced  than  at  a  mode  sparsity 
(a  poor  modal  density)  region  [17].  On  the  other  hand,  when  the  modal  overlap 
parameter  (Zy)  entertains  values  that  are  small  compared  with  imity,  the  levels  as 
a  function  of  (y),  tend  to  fluctuate.  The  fluctuations  are  pronounced  at  and  in  the 
vicinity  of  the  resonance  frequencies  of  those  satellite  oscillators  to  which  these 
small  values  of  (Zy)are  assigned.  At  and  in  the  vicinity  of  the  resonance 
frequencies  of  those  satellite  oscillators  to  which  (b^)  are  assigned  values  that 
approach  and  exceed  unity,  no  such  fluctuations  are  present;  e.g.,  see  Figs.  10c 
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and  11c,  and  Figs.  lOd  and  lid  and  contrast  them,  respectively  [18].  When 
variations  in  both  parameters  are  combined,  both  characteristics  can  be  identified 
in  the  levels  of  the  induced  loss  factor  e.g.,  see  Figs.  10  and  11  and 

contrast,  in  particular.  Figs.  lOe  and  1  le  with  Figs.  lOf  and  1  If. 
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Appendix  A 

The  sweeping  assumptions  rendered  in  Eq.  (15),  which  leads  to  Eqs.  (9d) 
and  (16),  culminating  in  Eq.  (18),  may  be  introduced  more  gradually.  The 
purpose  for  this  Appendix  is  to  effect  such  a  gradual  introduction.  In  this  manner 
when  some  of  these  assumptions  are  relieved,  reevaluation  of  the  induced  loss 
factor  rig{y)  may  be  readily  instituted. 

From  Eqs.  (9c)  and  (12),  the  design  demands  that 

{ccj.  +  a^r)  (x^)^  =  (l+m^^)  for  F=(l/2),  (Al) 

and  if  (x^)  becomes,  by  design,  unity  at  F=(l/2) ,  the  expression  reduces  to 

(a^  +«cr)  =  (1+^cr )  ^  =  (1/2) .  (A2a) 

From  Eq.  (12)  it  is  observed  with  satisfaction  that  (x®)  is,  indeed,  unity  at 
F  =  (1/2) .  On  the  other  hand,  from  Eqs.  (10)  and  (11)  one  obtains 

(x^f  TJr  =  {br  12)  [a{x^  {x°rf}ldr  ];  Tier  =  (b^r  /2)  [5{ac^(x°)^}/ar] ; 

jjr,  =[2iar+acr)ix^)^T^  {br[d{a,{x^rf}ldr]  +  bcr{d{acr{x^f}ldr]),  (A3) 

where  again,  (r)  is  allowed  to  have  a  continuous  connotation  as  explained  in 
Reference  1.  In  particular,  if  the  spring  factors  (a^)  and  (a^r)  are  independent 


of  (r) ,  Eq.  (A3)  simplifies  to  read 
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7,  =  (i,  ln)\riR)(x°f  ];  r,„  =  (b„  I  I 

Vrr  ~  Pr  Vcr  »  ~  ^  j  ^cr  ~ 

where  (x°)  is  extrapolated  and  interpolated  to  become  a  continuous  ftinction  of 

ir) 

[(ain(i;/5r)]  =  («  +  !)■'  lr(R)(x?fl  (A5) 

(A*<t)  =  («  +«cr'  [(*r«)+(*cr“c)]-  (A6) 

J,=[(;rA,)(ff +  !)->];  F,,  =  [(;r6,,)(fl  +  l)-‘ ],  (AT) 

and  from  Eq.  (A2a),  by  design 

(«+ac)  =  (l+^c)-  (A2b) 

From  Eqs.  (10)-(12)  and  (A2)-(A7)  one  may  cast  the  designed  two-vector  in  the 
form 

{(x„),  (r,„)}(.R)  =  {(l  +  m,)‘'^  (yff,)}  .  {(x”  ),  (n^r))  (R),  (AS) 

{(.X°„),  (11%))  (R)  =  {(x°r),  (bcr  !  !c)[y{R)(x°f  ]},  (A9) 

where 

[(l+^c)"’ [a(6^/^>cr)+«c]  ;  {AXQd) 


Pr  = 


=  P= 


(1  +  w^)  ^[a{blbc)+ac\,  b^=b,  bcr  =  bc  {AXOb) 


;  b  =  b^=b^^. 


1 


(.410c) 
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the  quantities  (w^),  (x®),  /(R)  and  (F^^)  are  stated  in  Eqs.  (Id),  (12),  (13)  and 

(A7),  respectively,  and  the  superscript  (o)  in  (x^)  and  (77^^)  recognizes  that 

(x%)  and  {tj%)  are  restricted  to  specific  impositions,  [cf.  Eq.  (4).]  For  the  sake 
of  computational  and  interpretive  advantage  the  validity  of  Eqs.  (A4)  and  (A  10c) 
is  universally  adopted  in  this  paper,  [cf  Eq.  (15).]  Under  this  imposition,  Eqs. 
(A8)  and  (A4)  simplify 

{(x„).(7„)}(/f)  =  {(l  +  m<,)''^(l?).7('0};  ft=[(^6)(J!  +  ir'],  (All) 

ri(7)  =  n,  =  ricr  =  ri„\  m  =  (A12) 

where  (x°)  is  stated  in  Eq.  (12),  and  becomes  equal  to  unity,  [cf  Eqs. 


(18a)  and  (16).] 
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Fig.  la.  A  master  oscillator  attached  to  a  set  of  sprung-masses. 


Fig.  lb.  A  set  of  satellite  oscillators  coupled  to  a  master  oscillator. 


Fig.  2.  (1.)  Resonance  frequency  distribution  (x^^)  and  (2.)  Corresponding  as¬ 
signed  loss  factor  (77^^ )  for  the  satellite  oscillators  as  a  function  of  the  normalized 
index  (F) ,  respectively.  (*) 

a.  1  and  a.2.  With  R  =  27  and  =  0.  [x^^  =  x"  ,  77^^  =  ri°^ .] 

(*)  When  regions  of  curves  clearly  overlap,  the  color  of  the  one  with  the  higher 
overlap  parameter  (b)  wins. 
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Fig.  2.  (1.)  Resonance  frequency  distribution  (x^^)  and  (2.)  Corresponding  as¬ 
signed  loss  factor  (77^^)  for  the  satellite  oscillators  as  a  fiinction  of  the  normalized 


index  (r ),  respectively.  (*) 

b.  1  and  b.2.  As  in  a.  1  and  a.2,  except  that  (R)  is  changed  from  (27)  to  (7). 
lx  =  x^ 


Fig.  2.  (1.)  Resonance  frequency  distribution  (x^^)  and  (2.)  Corresponding  as¬ 
signed  loss  factor  {-n,,)for  the  satellite  oscillators  as  a  function  of  the  normalized 
index  (F) ,  respectively.  (*) 

I 

c.l  and  C.2.  As  in  Fig.  2a,  except  that  (m^)  is  changed  from  (0)  to  (0.75). 


a. 


b. 


Fig.  3.  Induced  loss  factor  7],{y)  [=  ri°{y)],  as  a  function  of  (y),  for  a  stiffiiess 
control  coupling  form  with  =1.0  [«=0.0.],  g  =  =  0.  [Sprung-masses.] 

a.  With>  =  (0.1) 


b.  With  b  =  (2.0) 


Fig.  3.  Induced  loss  factor  ri^iy)  [=  7^  (y)],  as  a  function  of  (y),  for  a  stiffiiess 
control  coupling  form  with  =1.0  [a =0.0.],  g  =  =  0.  [Sprung-masses.] 

c.  With  ^  =  (10) 

d.  A  superimposition  of  b  =  (0.01),  (0.1)  and  (0.3),  and  (FOA).  (*) 
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e. 


f. 


Fig.  3.  Induced  loss  factor  rj^iy)  [=  7® (>-)],  as  a  function  of  (y),  for  a  stif&iess 
control  coupling  form  with  =1.0  [a: =0.0.],  g  =  m^=0.  [Sprung-masses.] 

e.  A  superimposition  of  a,  b,  and  c,  and  (FOA). 
f  As  in  3.e  above  except  that  R  =  (7)  and  not  (27). 


Fig.  4.  Induced  loss  factor  ii^iy),  as  a  function  of  (y),  for  a  stiffiiess  control 


coupling  form.  [R  =  27  and  (M^  /M)  =  0. 1 .] 


a.  -  0.75  \a  -  0.25.],  g  =  =  0.  [Strong  coupling.] 


b.  -  0.15  [a  —  0.85.],  g  =  =  0.  [Moderate  coupling.] 
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Fig.  6.  Induced  loss  factor  as  a  function  of  (y),  for  a  mass  control 

coupling  form.  [R  =  27  and  (M^  /M)  =  0. 1 .] 

~  S  —  =  1.15.],  =  0.15.  [Moderate  coupling.] 

^c~  S  —  =  1.03.],  —  0.03.  [Weak coupling.] 


Fig.  7.  Induced  loss  factor  as  a  function  of  (y),  for  mixed  control 

coupling  forms.  [R  =  27  and  (M^  /M)  =  0. 1 .] 


a.  a^.  -  -  0.75  [a  -  1.0.],  g  =  0 .  [Very  strong  coupling.] 


b.  -  0.53  [a  —  0.47.],  g  =  0.54,  =  0 .  [Strong  coupling.] 


Fig.  7.  Induced  loss  factor  T]^{y),  as  a  function  of  (y),  for  mixed  control 
coupling  forms.  [R  =  27  and  {M^  I M)  =  0.1.] 

c.  -  0.10  [a  =  0.90.],  g  =  0.11,  [Moderate  coupling.] 

d.  -  0.02  [a  =  0.98.],  g  =  0.022,  =  0.  [Weak  coupling.] 


Fig.  8.  Pseudo-statistical  variations: 


Fig.  9.  Modification  of  pseudo-statistical  variations: 
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Fig.  10.  Induced  loss  factor  7,(3;)  [=7f(:);)],  as  a  function  of  (:);),  for  a  strong  coupling  defined  by  =  1.0 [a  =  0.0.], 
=  f  =  0;  implementing  a  number  of  pseudo-statistical  variations.  [R  =  27  and  /M)  =  0.1 .] 

First  order  approximation  is  superimposed. 


Fig.  11.  As  in  Fig.  10,  except  that  the  coupling  is  a  mix  of  stif&iess  and  gyroscopic  control:  =  0.1[o'=0.9.], 

g  =  0.1 1,  =  0.  The  coupling  strength  here  is  moderate. 


Fig.  11.  As  in  Fig.  10,  except  that  the  coupling  is  a  mix  of  stiffiiess  and  gyroscopic  control:  =  0.1[a=0.9.], 

g  =  0.1 1,  =  0.  The  coupling  strength  here  is  moderate. 
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